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ABSTRACT 

qq ■ We investigate formation processes and physical properties of globular cluster 

00 | systems (GCSs) in galaxies based on high-resolution cosmological simulations with 

, globular clusters. We focus on metal-poor clusters (MPCs) and correlations with their 

£S| ' host galaxies by assuming that MPC formation is truncated at a high redshift (ztrun > 

| 6). We find that the correlation between mean metallicities (Z gc ) of MPCs and their 

■ host galaxy luminosities (L) flattens from z — z trun to z = 0. We also find that 

the observed relation (Z gc oc L 015 ) in MPCs can be reproduced well in the models 

. with Z gc ~ L 5 at z = z tIU n when z trun ~ 10, if mass-to-light-ratios are assumed 

to be constant at z — z tlun . A flatter L — Z gc at z = z trun is found to be required 
to explain the observed relation for constant mass-to-light-ratio models with lower 
z = ztrun- However, better agreement with the observed relation is found for models 
with different mass-to-light-ratios between z = Ztrun and z = 0. It is also found that 
the observed color- magnitude relation of luminous MPCs (i.e., "blue tilts") may only 
r* ■ have a small contribution from the stripped stellar nuclei of dwarf galaxies, which have 

. ^ ' nuclei masses that correlate with their total mass at z = ztrun- The simulated blue 

, tilts are found to be seen more clearly in more massive galaxies, which reflects the 

■ fact that more massive galaxies at z — are formed from a larger number of dwarfs 

with stellar nuclei formed at z > z trun . The half-number radii (R c ) of GCSs, velocity 
dispersions of GCSs (er), and their host galaxy masses (M n ) are found to be correlated 
with one another such that R c oc Mh°' 57 and a cx M n ' 32 . Based on these results, we 
discuss the link between hierarchical merging histories of galaxies and the physical 
properties of MPCs, the origin of the L — Z gc relation, and non-homology of GCSs. 

Key words: globular clusters: general - galaxies: star clusters - galaxies: evolution - 
galaxies: stellar content 
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1 INTRODUCTION 

A growing number of observational studies of globular clus- 
ter systems (GCSs) have revealed interesting correlations 
between physical properties of GCSs and those of their host 
galaxies (see Brodie & Strader 2006 for a recent review). 
For example, Strader et al. (2004) found that mean colors 
of metal-poor globular clusters (GCs) correlate with lumi- 
nosities of their host galaxies and accordingly suggested that 
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mean GC metallicities (Z gc ) depend on luminosities (L) of 
their host galaxies such that Z sc cx £° 15 ± 03 Similar re- 
sults (Zg C oc £°- 16 ± ' 04 ) were found using the Advanced 
Camera for Surveys on the Hubble Space Telescope (HST) 
by Peng et al. (2006; P06). Brodie & Strader (2006) have 
suggested that the relation was steeper in the past. They 
argue that more enriched metal-poor GCs form first in low 
mass objects (from early collapsing peaks at high redshift). 
These low-mass building blocks merge, ultimately forming 
the massive galaxies of high density regions today. Thus the 
slope, particularly at the high mass end, flattens over time. 

Although such observed correlations between the physi- 
cal properties of GCSs and their host galaxies have been sug- 
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gested to contain fossil information of galaxy formation and 
evolution (e.g., Harris 1991; Forbes & Forte 2001; West et 
al. 2004; Brodie & Strader 2006), the details remain largely 
unclear owing to the lack of theoretical and numerical stud- 
ies of GCSs. Based on dissipationless numerical simulations 
of major galaxy mergers with GCs, Bekki & Forbes (2006) 
first demonstrated that the GCS radial density profiles de- 
pendent on galaxy luminosity can be understood in terms 
of the number of major merger events experienced. The de- 
pendence of GC destruction on host galaxy mass has also 
been suggested to be important for understanding the in- 
ner radial density profiles (Baumgardt 1998; Vesperini et al. 
2003; Capuzzo-Dolcetta 2006). 

The formation of GCs in low-mass dark matter halos 
at high redshifts have been investigated by numerical and 
theoretical studies (e.g., Broom & Clarke 2002; Mashchenko 
& Sills 2005). Furthermore, the physical properties of GCSs 
have recently been discussed by several authors based on hi- 
erarchical galaxy formation scenarios where GC formation 
occurs at high redshifts (e.g., Beasley et al. 2002; Kravtsov 
& Gnedin 2005; Santos 2003; Bekki 2005; Rhode et al. 2005; 
Yahagi & Bekki 2005; Moore et al. 2006; Bekki et al. 2006; 
Bekki & Yahagi 2006). For example, Beasley et al. (2002) 
first demonstrated that the observed bimodal color distribu- 
tions of GCSs in elliptical galaxies can be reproduced by a 
semi-analytic galaxy formation model. They also presented 
GC colour vs galaxy magnitude relations, however they did 
not fully reproduce the observed trend. Although these pre- 
vious works discussed some physical properties of GCSs in 
galaxies, they did not explore in detail. The correlations be- 
tween GCSs and those of their host galaxy. Thus it remains 
unclear how the physical properties of GCSs evolve in a hi- 
erarchical merging cosmology. 

The properties of metal-poor GCs (MPCs) ([Fe/H] ~ 
— 1.5) such as their very old ages, low metallicities and ex- 
tended spatial distribution all suggest formation at early 
times, when the low mass building blocks of galaxies formed. 
The purpose of this current paper is to investigate the physi- 
cal properties of MPCs and their scaling relations with their 
host galaxies based on high-resolution cosmological simu- 
lations that follow realistic merging and accretion histo- 
ries of galaxies. We investigate the structural, kinematical, 
and chemical properties of MPCs. We also present several 
predictions for the expected correlations between proper- 
ties of MPC systems (such as half-number radii, effective 
surface number densities, and velocity dispersions). Since 
the present model is based on dissipationless simulations of 
MPC formation within halos at high redshifts, the origin 
of metal-rich GCs (MRCs) ([Fe/H] ~ -0.5) formed some- 
what later during dissipative merging will not be discussed. 
We plan to investigate MRCs in a future paper by combin- 
ing the present high-resolution cosmological simulation with 
semi-analytic models similar to those used by Beasley et al. 
(2002). 

The plan of this paper is as follows: In the next sec- 
tion, we describe our numerical models of MPC formation 
and their assumed initial properties. In §3, we present our 
numerical results on (i) L — Z sc relations, (ii) scaling re- 
lations between properties of MPCs, and (iii) metallicity- 
magnitude relations for luminous MPCs in massive galaxies 
(the so-called 'blue-tilt'). In §4, we discuss our results. We 
summarize our conclusions in §5. Appendices explore scaling 



relation dependency on Ms/L and z irun . Throughout this 
paper, GCSs in our simulations are composed only of MPCs. 



2 THE MODEL 

2.1 Identification of MPCs 

We simulate the large scale structure of GCs in a ACDM 
Universe with Q. = 0.3, A = 0.7, H = 70 km s" 1 Mpc" 1 , 
and a s = 0.9 by using the Adaptive Mesh Refinement 
N— body code developed by Yahagi (2005) and Yahagi et 
al. (2004), which is a vectorized and parallelized version 
of the code described in Yahagi & Yoshii (2001). We use 
512 3 collisionless dark matter (DM) particles in a simula- 
tion with the box size of 70/i _1 Mpc and the total mass of 
4.08 x 10 16 M Q . We start simulations at z = 41 and follow 
it until z = in order to investigate physical properties of 
old GCs outside and inside virialized dark matter halos. We 
used COSMICS (Cosmological Initial Conditions and Mi- 
crowave Anisotropy Codes), which is a package of fortran 
programs for generating Gaussian random initial conditions 
for nonlinear structure formation simulations (Bertschinger 
1995, 2001). 

Our method of investigating GC properties is described 
as follows. Firstly, we select virialized dark matter subha- 
los at a truncation redshift z — zt run by using the friends- 
of-friends (FoF) algorithm (Davis et al. 1985) with a fixed 
linking length of 0.2 times the mean DM particle separation. 
The minimum particle number iV m i n for halos is set to be 
10. For each individual virialized subhalo with a half-mass 
radius of i? n , particles within i? n /3 are labeled as "GC" par- 
ticles and are considered to be old MPCs. This procedure 
for defining GC particles is based on the assumption that 
energy dissipation via radiative cooling allows baryons to 
fall into the deepest potential well of dark-matter halos and 
finally to be converted into GCs. The value of the trunca- 
tion radius (-Rtr.gc = Rh/3) is chosen, because the size of the 
old GCs in the Galactic GC system (i.e., the radius within 
which most Galactic old GCs are located) is similar to -Rh/3 
of the dark matter halo in a dynamical model of the Galaxy 
(Bekki et al. 2005). We assume that old MPC formation is 
truncated completely after z = ztrun and investigate the de- 
pendences of the results on zt run . Physical motivation for 
the truncation of GC formation is described later. 

Secondly, we follow GC particles formed before z = 
ztrun until z = and thereby derive locations (x,y,z) and 
velocities (v x ,v y ,v z ) of GCs at z = 0. We then identify viri- 
alized halos at z = with the FoF algorithm and investigate 
whether each GC is within a halo. If GCs are found to be 
within a halo, the mass of the host halo (M n ) and physical 
properties of the GCS are investigated. If a GC is not in 
any halo, it is regarded as an intergalactic GC. The num- 
ber fraction of these depends on ztrun but is typically less 
than 1 % (e.g., ~ 0.3 % for = 10). We don't consider 
intergalactic GCs further in this current paper. 

Thus, the present simulations enable us to investigate 
the physical properties only for old MPCs due to the adopted 
assumption of collisionless simulations. The physical proper- 
ties of MRCs which may form later during dissipative merger 
events (e.g., Ashman & Zepf 1992) are not investigated. 
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Figure 1. Spatial distributions of all MPCs projected onto the x-y plane in the large-scale cosmological simulation (left) and spatial 
distributions of MPCs within the selected four galaxy-scale halos with similar masses (right) for the model with ztrun = 10. The halo 
masses M h are 9.3 X 1O 12 M , 9.6 X 1O 12 M , 6.1 X 10 12 M Q , and 8.1 X 1O 12 M , for the upper left (a), the upper right (b), the lower left 
(c), and lower right (d), respectively. The four circles shown by solid lines represent the locations of the four halos in the left panel. The 
horizontal bar is 20 kpc. The total number of MPCs range from 31 to 249 in the four halos, which reflects the different merging histories 
with MPCs formed before ztrun =10. 



2.2 Initial properties of MPCs in galaxy-scale 
halos 

2.2.1 Radial density profiles 

We assume that the initial radial profiles of GCSs (p(r)) in 
subhalos at z = ztrun are the same as those described by the 
universal "NFW" profiles (Navarro, Frenk & White 1996) 
with p(r) oc r~ 3 in their outer parts. The mean mass of sub- 
halos at z — zt rU n in the present simulations is roughly sim- 
ilar to the total mass of dwarf galaxies today. Minniti et al. 
(1996) found that the projected (R) density profiles of GCSs 
in dwarfs is approximated as p(R) oc R~ 2 , which translates 
roughly to p(r) cx r~ 3 using a canonical conversion formula 
from p(R) into p(r) (Binney & Tremaine 1987). Therefore, 
the above r dependency can be regarded as reasonable. 
Below we mainly show the fiducial model with p(r) similar 
to the NFW profiles and i? tr , gc = R h /3. 

Although we base our GC models on observational re- 
sults of GCSs at z=0, we can not confirm whether the above 
p(r) and Ru tgc values are reasonable for GCSs in low-mass 
subhalos at z > 6 owing to the lack of observational data 
for GCSs at high redshifts. Although our present results on 
the L — Z gc relation and the blue tilt at z — do not de- 
pend strongly on the initial distributions of MPCs, struc- 
tural properties of GCSs do depend on the spatial distribu- 
tion. For example, if we choose smaller Ru, gc at z — z trU n, 
the final projected number distributions of GCSs becomes 
more compact. 



2.2.2 The L-Z sc relation 

In order to investigate the mean metallicities of MPCs in 
galaxy-scale halos at z = 0, we need to allocate initial metal- 
licities to all GCs formed before z = Ztrun- We assume that 
GCs in a halo at z trU n have identical metallicities of Z gc and 
that the metallicities are a function of the total stellar mass 
(or luminosity) of the halo. Lotz et al. (2004) have found that 
Z gc oc L a&c , where a gc = 0.2 for B— band luminosities in dEs 
and dE,Ns. This relation is similar to Z gc cx L discovered 
by Strader et al. (2004). Dekel & Silk (1980) demonstrated 
that the stellar metallicities (Z) of dwarf galaxies embed- 
ded in massive dark matter halos correlate with luminosi- 
ties (L) of the dwarfs such that L oc Z 2 ' 7 . This translates 
to Z gc tx L ' 37 if GCs follow field stars in dwarfs. Prompted 
by these theoretical and observational studies, we adopt the 
following relation at high redshift between [Fe/H] (or log 
Z gc ) of MPCs and L: 

[Fc/H]=a gc xlogL + /3 gc , (1) 

where a gc and f3 gc are set to be 0.5 and -6.0 in the present 
study. We show the models with these values, because they 
are more consistent with the observed L — Z gc relation. 

In order to derive L for halos with masses M n at 
z = 0, we assume a mass-to-light-ratio (M h /L) dependent 
on Mh. Recent observational studies based on galaxy lumi- 
nosity functions for luminous galaxies suggests that Mh/L 
depends on M h as M h /L oc M h 33 (Marinoni & Hudson 
2002). Dekel & Silk (1986) proposed that M h /L cc L~ 0,37 
for low-luminosity dwarf galaxies, which can be interpreted 
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as Mh/L oc M h ~ 0,55 . These studies suggest that the -/In- 
dependences of Mh/L are different between low- and high- 
luminosity galaxies (see also Zaritsky et al. 2006). We thus 
adopt two different Mh-dependent Mh/L ratios for galaxy- 
scale halos in the present study. For halos above a threshold 
halo mass of M h , t h, we adopt the following: 

M h /L = Cml x {^)'\ (2) 

where Cml is a constant and the value of Mh/L at Mh = 
M h ,th- For halos below a threshold halo mass of M h ,th, we 
adopt the following: 

Mh/L = Cml x ( — , (3) 

JWh.th 

We assume that Cml = 10 and M h ,th = 1O 11 M are rea- 
sonable values. Since the threshold mass M n ,th has not yet 
been so precisely determined by observational studies (e.g., 
Zaritsky et al. 2006), M h ,th can be a free parameter in our 
simulations. We however confirm that this Mh,th for a rea- 
sonable range is not so important as other parameters (e.g., 
Mh/L dependences). We thus mainly show the results of the 
models with M h ,th = 10 n M Q . 

In calculating Z gc from L at z = z trU n, we need to de- 
rive L from Mh at z = z trU n- It is, however, observationally 
unclear what a reasonable Mh-dependent Mh/L is for ha- 
los at z = 2trun- We accordingly investigate two extreme 
cases: One is that Mh/L = Cml for all halos at z = Ztrun 
and the other is that the Mh-dependence of Mh/L is the 
same between z = (equations 2 and 3) and z — ztrun- The 
total stellar mass (M a ) of a halo is estimated from Mh by 
using the adopted Mh-dependent Mh/L and the following 
stellar-mass-to-light- ratio (M B /L). 

It is highly likely that M h -dependences of M h /L are 
different between different redshifts, because of (i) evolution 
(e.g., aging) of stellar populations with different ages and 
metallicities and (ii) evolution of baryonic mass fraction. 
Although observational studies on the redshift evolution of 
the "Fundamental Plane" of early type galaxies (e.g., van de 
Ven et al. 2003) have provided some clues to stellar popu- 
lation evolution it remains observationally unclear whether 
evolution of baryonic mass fraction is seen in galaxies: kine- 
matical data sets for precisely estimating total masses of 
galaxies are currently unavailable for very high- 2 galaxies. 
The present models, which do not allow us to make explicit 
predictions on redshift evolution of Mh/L, have difficulties 
to determine which of the above two is responsible for the 
possible redshift evolution in M h -dependences of Mh/L. We 
therefore do not intend to discuss the origin of the redshift 
evolution in Mh-dependences of Mh/L in the present paper. 

The stellar population synthesis models by Vazdekis et 
al. (1996) for a Salpeter IMF show that M s /L in the K-band 
for stellar populations with [Fe/H] = -1.7 (0.0) is 0.22 (0.39) 
for 0.5 Gyr and 2.03 (4.84) for 12.6 Gyr. Although M s /L 
depends strongly on age and metallicity of the stellar pop- 
ulation, we use a constant M s /L of 1 for all halos because 
of the lack of information on stellar populations in the sim- 
ulated galaxy-scale halos. A discussion on the dependences 
on M s /L is given in Appendix B. 

Given the fact that most of the virialized halos at ztrun 
are of low- mass (< 1O 1O M0) then the low-metallicities used 
to derive M B /L should be reasonable. We adopt M s /L = 1, 



which is the average value between 0.5 and 12.6 Gyrs in 
metal-poor stellar populations both for Ztrun and z — for 
halos with different masses. Although this adoption is more 
idealized, it helps us to derive more clearly the evolution 
of the M s — Z S c relation between z tIun and 2 = with- 
out having complicated dependences of M B /L on ages and 
metallicities of stars in halos (that can not be modeled di- 
rectly in the present study). If we adopt a higher M B /L (e.g., 
~ 4), the simulated M s — Z gc relation at z = is shifted to 
lower metallicities for a given stellar mass but with the same 
relation slope. 



2.2.3 GC Luminosity function 

The luminosity function of GCs has been observationally 
suggested to be universal and have the following form (e.g., 
Harris 1991): 

$(M) = Co x exp(-(M - M£) 2 /2a m 2 ), (4) 

where Co is a constant, My = —7.23 ±0.23, and er m = 1.25 
mag (Harris 1991). We adopt this luminosity function and 
allocate luminosities to GC particles by generating random 
numbers for —11 < My < —4. We note that our results are 
not very sensitive to the choice of M v or cr m . 



2.3 Initial properties of galactic nuclei 

Previous studies have suggested that stellar galactic nu- 
clei (GN) in nucleated low-mass galaxies might be observed 
as very massive GCs after the host galaxy of the GN is 
destroyed during merging with, and accretion onto, more 
massive galaxies (e.g., Zinnecker et al. 1988; Freeman 1993; 
Bassino et al. 1994; Bekki & Freeman 2003; Bekki & Chiba 
2004). The physical properties of very massive GCs in the 
Galaxy and M31 (e.g., lu Cen and Gl) are observed to be 
different from those of "normal" GCs (e.g., Freeman 1993). 
We identify stellar galactic nuclei particles (referred to as 
GN particles) at z = z trU n and follow their evolution until 
z = 0. 

The ACS Virgo Cluster Survey by Cote et al. (2006) 
recently found that the mass fraction of GN (stellar galactic 
nuclei) to their host galaxies is typically 0.3% in early-type 
galaxies. We thus adopt the following relation between the 
total masses of GN (M gn ) and the total stellar masses of 
their hosts M s : 

M gn = / gn X M s , (5) 

where / gn is considered to be a free parameter in the present 
study. Recent observations of stellar populations of SGN 
in nearby galaxies suggest that GN have a significant con- 
tribution from young stellar populations and thus suggest 
that GN are growing with time (e.g., Walcher et al. 2006). 
Therefore it is reasonable to adopt an / gn value significantly 
smaller than today's 0.3%. We generally adopt / gn = 0.001 
but also discuss the dependence of our present results on 

/gn. 

Lotz et al. (2004) and Cote et al. (2006) found a corre- 
lation between M gn (or GN luminosities) and colors of GN, 
which suggests a mass-metallicity relation for GN (i.e., more 
massive GN are more metal-rich). This correlation and the 



above M gn -M s relation implies that more luminous galax- 
ies have more metal-rich GN. Lotz et al. (2004) found that 
the metallicities of field stars of dE,Ns (Zs) correlates with 
galaxy luminosities (L) such that Zs oc L 0,4 . Kravtsov & 
Gnedin (2005) showed that galaxies in their cosmological 
simulations exhibit a strong correlation between the stel- 
lar mass (Ms) and the average metallicity (Z) of stars (de- 
scribed as Zs oc Ms ' 5 ). Guided by these studies, we assume 
the following relation between GN metallicities (log Z gn ) 
and luminosities of their hosts (L) : 

[Fe/H] = a gn x log L + /3 gn , (6) 

where a gn and /3 gn are set to be the same as those adopted 
for MPCs denoted as a gc and /3 gc (i.e., 0.5 and -6.0, re- 
spectively) in the present study for consistency. We discuss 
whether the adopted models can provide a physical explana- 
tions for the origin of mass-metallicity relation of luminous 
GCs around massive galaxies in §4.2. 

2.4 Truncation of MPC formation 

Semi-analytic galaxy formation models based on a hierar- 
chical clustering scenario have shown that models with a 
truncation of GC formation at z ~ 5 can better repro- 
duce the observed color bimodality of GCs in early-type 
galaxies (Beasley et al. 2002). Recent cosmological simula- 
tions have demonstrated that truncation of GC formation 
by some physical mechanism (e.g., reionization) is necessary 
to explain the very high specific frequency (Sn) in cluster 
ellipticals, structural properties of the Galactic old stars and 
GCs, and the mass-dependent Sn trend (Santos 2003; Bekki 
2005; Bekki & Chiba 2005; Rhode et al. 2005; Moore et al. 
2006; Bekki et al. 2006). Observational studies however have 
not yet found any strong evidence for the truncation of GC 
formation in GCSs of nearby galaxies. 

If •Ztrun is closely associated with the completion of cos- 
mic reionization, then z tr un may wen range from 6 (Fan et 
al. 2003) to 20 (Kogut et al. 2003). Although we investi- 
gate models with different z tIU n ( i.e. 6, 10, and 15), here 
we present the results of the model with Ztfun = 10. This 
is because the latest observations by Microwave Anisotropy 
Probe ( WMAP) suggested that the epoch of reionization 
is 2 = 10.9^2*3 (Page et al. 2006), and because numerical 
simulations suggest that models with higher Zt run (> 10) 
are more consistent with observations of MPCs (e.g., Bekki 
2005; Moore et al. 2006). 
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Figure 2. Mean metallicities ([Fe/H]) of metal-poor GCSs in 
simulated galaxies as a function of the total halo masses of the 
galaxies (log^QM n ) for the model with 2trun 

= 10. Each small 

dot represents a galaxy with a GCS at z = 0. The thin blue line 
represents the initial relation (z = 10) between Z gc and M h , where 
Z gc is the mean metallicity of their MPCs. The thin green line 
represents the mean value of the simulated relation today (z=0). 
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2.5 Galaxy-scale halos 

We simulate the structural, kinematical, and chemical prop- 
erties of 114596 GCs formed in 25096 low- mass halos before 
ztrun = 10. The physical properties of GCSs in 2830 halos 
with masses 4.0 x 1O 9 M < M h < 6.5 x 10 14 M Q at z = 
are analyzed and their correlations with host properties are 
investigated. Since we are interested in the GCS proper- 
ties in galaxy-scale halos we only show the results for GCSs 
in halos with masses lO n M < M h < 1O 13 M (and with 
1 < N gc < 246). 

We investigate the relations between three properties of 
GCSs with iV gc > 4, i.e. the half-number radii of projected 
GC distribution (Re), the surface number densities at R e 



_2 Uj i I i i i i I i i i i I i i i i I i i i lJ 

9 10 11 12 13 

l°g,0 M s ( M Q ) 

Figure 3. Mean metallicities ([Fe/H]) of metal-poor GCSs in 
simulated galaxies as a function of the total stellar masses of the 
galaxies (log 10 Ms) for the model with ztrun = 10. Each small 
dot represents a galaxy with a GCS at z = 0. Here we use a 
variable M/L at z = and a constant M/L at z = 10. The thin 
blue line represents the initial Z gc oc L° relation, where L is 
the total luminosity of galaxies and Z gc is the mean metallicity 
of their MPCs. The thin green line represents the mean value of 
the simulated relation today (z=0). For comparison, the observed 
L — Z gc relation of Z gc oc L 01B (P06) is shown by a thick red 
line. 



6 K. Bekki, H. Yahagi, and D. A. Forbes 



CD 

a 
p 



1.2 



0.8 



CD 
N 



0.6 



0.4 



0.2 



-0.5 



-1.8 



-1.6 



-1.4 



-1.2 



[Fe/H] 



Figure 4. The metallicity distribution function (MDF) of GCSs 
from progenitor low-mass halos at z = 10 that forms a galaxy- 
scale halo of M h = 5.0 x 1O 12 M at z = 0. This halo with a 
mean Z gc = —1.39 (dashed line) at z = originates from a low- 

-1.20. 
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The GCS at 2 = forms via hierarchical merging and accretion, 
between z = 10 and z = 0, of 74 low-mass halos each with a GCSs 
containing a range of metallicities. 
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Figure 6. The same as Figure 3 but for the model with 2t rU n 
10 and variable M/L both at z = and at z = 10. 
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Figure 5. The metallicity difference between GCSs in galaxy- 
scale halos at 2 = and those in their progenitor halos at z = 10 
(<5Z gc = Z gc (z = 10) — Z gc (z = 0)) with halo masses (Mh) at 
2 = 0. Small dots represent GCSs in galaxy-scale halos and big 
filled circles arc the mean values of SZ gc . 



(J ), and the line-of-sight velocity dispersion for all MPCs 
in each halo (a). We derive these three properties, because 
they are observationally feasible to measure. Since M s /L 
is fixed at a constant value of 1 in the present study, the 
simulated M s — Z gc relation is virtually the same as the 
L — Z gc relation that is often discussed. The dependences of 
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3 RESULT 

3.1 The L - Z gc relation 

Fig. 1 shows the large-scale (~ 50 Mpc) distribution of 
MPCs at z — and MPCs within galaxy-scale halos with 
similar halo masses (Mh ~ 8 x 1O 12 M0) in the model 



with Ztrun = 10. It it clear from Fig. 1 that the spatial 
distribution of MPCs in galaxy-scale halos are quite di- 
verse. For example, the GCS distribution in the halo with 
Mh = 6.1 x 1O 12 M0 appears to be more flattened than that 
in the halo with M h = 9.6 x 1O 12 M . The mean spatial 
density of MPCs within the central 50 kpc of the halo with 
M h = 8.1 x 1O 12 M is a factor of ~ 8 higher than that 
of the halo with M h = 9.3 x 10 12 M Q . These diverse GCS 
structures are due to the differences of merging histories 
of halos (e.g., the number of low-mass subhalos that were 
virialized before z tV un which merged to formed a halo). The 
mean metallicities and metallicity distribution functions of 
GCSs in galaxies are also diverse depending on their merging 
histories. 

Fig. 2 shows how the initial Mh — Z gc relation evolves 
between z = 10 and z — for the model with Ztmn = 10. 
Although the dispersion in [Fe/H] at z = is significantly 
large for a given metallicity bin, the Afh — Z gc relation at 
z = 10 becomes flatter at z = 0. This result clearly indicates 
that the mean metallicities of MPCs in galaxies correlate 
with the total masses of their host galaxies and thus with 
those of dark matter halos. This also demonstrates that the 
simulated L — Z gc (discussed later) results from the Mh — 
Z gc relation rather than the adopted assumptions on Mh/L. 
The dependences of this correlation on ^trun are given in 
Appendix A. 

Fig. 3 shows that the initially steep L — Z gc relation 
at z = 10 becomes significantly flatter at z — due to 
hierarchical merging of low-mass galaxies formed at high z. 
Fig. 3 also shows that the dispersion in [Fe/H] for a given 
mass range of M s is quite large. Fig. 3 demonstrates that the 
observed L — Z gc relation (i.e., Z gc <x L 016 ) at z = can 
be reasonably well reproduced by the model with an initial 
L — Z gc relation of Z gc cx L ' 5 at z — 10. We find that such 
a steep initial L — Z gc relation at ztrun is indispensable for 
reproducing the very flat yet significant correlation between 
0. The flattening of the initial L — Z gc 
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relation can be seen also in models with different 2 trU n (see 
Appendix C) 

Figs. 4 and 5 show why the initial steep L — Z gc rela- 
tion becomes flatter during the hierarchical growth of galax- 
ies. Fig. 4 shows that the final (2 = 0) mean metallic- 
ity of the MPCs in a galaxy with M h = 5.0 x 1O 12 M 
is lower (Z gc = —1.39) than the initial mean metallic- 
ity (Z gc = —1.20) of its progenitor low-mass halo with 
Mh = 4.0 x 10 10 Mq, which is the most massive halo among 
those forming this galaxy at z = 10. This galaxy grows via 
merging and accretion of 74 low-mass halos, all with GCSs 
of different metallicities between z = 10 and 2 = 0. Since the 
MPCs in the accreted "building blocks" are metal-poor, the 
metallicity of the GCS in the final galaxy is lower than that 
of the GCS of its progenitor halo. There is thus a decrease 
in the mean metallicities of GCSs during galaxy evolution 
between 2t rU n and 2 = which is driven by merging and ac- 
cretion of low-mass halos and their lower metallicity GCSs. 

Fig. 5 shows the differences in mean metallicities (Z gc ) 
between GCSs in galaxies at 2 = and those in their pro- 
genitor low-mass halos at 2 = 10 (i.e. 8Z gc = Z gc (z = 
10) — Zgc (2 = 0)). Although the dispersion in SZ gc is quite 
large in more massive galaxies, the absolute values of 8Z gc 
are significantly larger in more massive galaxies. This means 
that more massive galaxies at 2 = have experienced a 
larger number of merging and accretion events of lower- 
mass halos with lower Z gc , so that the final Z gc at 2 = 
is lower. Low-mass galaxies at 2 = 0, on the other hand, 
have not experienced so many mergers and accretion events 
that act to increase their masses and decrease their Z gc . 
Such galaxies are of low-mass because they have only grown 
slightly through merging/accretion since their virialization 
at 2 = 10. As a result of this and forming relatively many 
GCs in the progenitor halo, their Z gc does not change so 
much from their original values at 2 = 10. Thus the mass- 
dependence of 5Z gc is the main reason for the flattening of 
the L — Zgc relation between 2 trU n and 2 = 0. 

Fig. 6 shows that if the M/L ratios in galaxies are vari- 
able at 2 = 10 and 2 = 0, then the simulated Z gc oc L 0,3 
relation is not a good fit to the observed one. These results, 
combined with those in Fig. 3, demonstrate that the model 
with a constant M/L at 2 = 10 can better reproduce the 
observed flat L — Z gc relation for a given initial L — Z gc rela- 
tion. Given the fact that both the M h -dependence of M/L 
and L — Z gc relation at 2 = 10 are observationally unclear, 
there are two possible interpretations for the derived steeper 
L — Z gc relation at 2 = in the variable M/L models. One 
possibility is that the variable M/L is more reasonable, and 
thus the initial L — Z gc relation at 2 = 10 should be much 
flatter than the adopted one of Z gc oc L . The other pos- 
sibility is that a constant M/L is more reasonable and the 
observed M h -dependence of M/L at 2 = is achieved during 
galaxy evolution between 2 = 10 and 2 = 0. 

An initial L — Z gc relation at 2 = 10 is required to have 
Zgc oc L ' to explain the observed relation in a variable 
M/L model, because merging flattens the slope by a factor 
of 1.7 in this model. This L — Z gc relation at 2 = 10 is flatter 
than the luminosity-metallicity relation of low-mass galaxies 
(Z oc L ' 37 ) predicted by Dekel & Silk (1980). Given the 
fact that there are no observational constraints on an initial 



constant M/L models with initially steep L — Z gc relations 
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Figure 7. Velocity dispersions (a, top), effective surface num- 
ber densities (I c , middle), and half-number radii (R e , bottom) 
of MPCs vs the total halo masses of galaxies at z = for the 
model with ztrun = 10. The solid line in each panel represents 
the least-square fit to the simulation data. 



are better than variable M/L models with less steep L — Z gc 
relations at this stage. 

The present simulations thus show two extreme cases 
of L — Zgc evolution depending on the assumed M/L at 
z — 2 t run • The present simulations however confirm that the 
initial L — Z gc relation of GCSs at 2 = 2 trU n becomes sig- 
nificantly flatter by 2 = 0, irrespective of different (yet rea- 
sonable) initial Mh-dependences of M/L and different 2 trU n. 
These results thus imply that the L — Z gc relation evolves 
with 2 in such a way that the L — Z gc relation becomes flatter 
at lower 2. 



3.2 Scaling relations 

Fig. 7 shows the scaling relations between three projected 
properties of GCSs, R e (half- number radii), 7 e (surface num- 
ber densities at Re), and a (line-of-sight velocity dispersion). 
Fig. 7 demonstrates that a and R c strongly correlate with 
M h in such a way that a and R e are both larger for larger 
Mh. I c however correlates weakly with Mh such that J c is 
smaller for larger Mh . The least square fits to the simulation 
data show that Re oc M h ' 57 , I e oc M h ~ - 65 , and a cc M h - 32 . 
These scaling relations appear to be different from those de- 
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Figure 8. The same as Figure 7 but for the dependences 
total stellar masses of galaxies. 



rived for dark matter halos (Kormendy & Freeman 2004), 
which suggests that GCSs in our models do not simply trace 
the distributions of dark matter halos. These scaling rela- 
tions reflect the fact that MPCs in galaxy-scale halos at 
z — originate from the central regions of low-mass halos 
at z = 10 and thus have different structures from those of 
dark matter halos at z = 0. 

Fig. 8 shows how the properties of GCSs (R c , J c a) 
correlate with total stellar masses (M a ) in halos and confirm 
that these correlations are similar to those derived in Fig. 
7. The least square fits to the simulation data show that 
R c oc M s - 85 , / c oc M B -°- 97 , and a oc M B 0AS . The derived 
M s -dependences are steeper than the Mh-dependences due 
to the adopted Mh variation with M/L in galaxies at z = 
0. These results suggest that future observational studies 
of GCSs for a statistically significant number of GCSs can 
assess the validity of the present MPC formation model. The 
dependency of the CCS scaling relations on Ztmn are given 
and their physical meaning discussed in Appendix D. 

Fig. 9 shows that the properties of GCSs (i.e., R c , 7 C , 
and a) in our models do not follow a relation for dynamical 
systems in virial equilibrium (i.e., R c = Cvpie _1 c 2 ). Given 
the fact that the simulated GCSs at z — are in dynam- 
ically equilibrium, this apparent deviation from the virial 
relation suggests structural and kinematical non-homology 
in GCSs, as observationally and theoretically suggested for 
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Figure 9. Half-number radii for projected MPCs as a function 
of the combination of I c (surface-number density at R e ) and a 
(line-of-sight velocity dispersion for MPCs) for the model with 
ztrun = 10. The solid line represents the "virial plane" defined as 
Rc = Cvp^c~ 1 o' 2 . Note that the distribution of GCSs deviates 
significantly from the virial plane. Physical interpretations of this 
result is given in the main text. 
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Figure 10. The same as Figure 9 but for the different abscissa. 
The solid line represents the "Fundamental Plane" defined as 
ij c = Cfp-Z o ~ ' 46 <t ' 60 . The dispersion of the simulation data 
along this line is small, suggesting that GCSs are located on a 
Fundamental Plane. Physical interpretations of this result is given 
in the main text. 
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Figure 11. Distribution of MPCs in the M v -[Fe/H] plane for 
the halo with M h = 3.0 X 1O 13 M and jV gn = 62 in the model 
with / gn = 0.001. Big dots represent the mean My for each of the 
five metallicity bins. This halo experienced merging/ accretion of 
nucleated galaxies 62 times and thus now has 62 stripped nuclei 
in its GCS. Note that a Mv-[Fe/H] relation can be clearly seen for 
MPCs with —12 < My < —9 (mag). This relation originates from 
the L — Z gn relation of low-mass "building blocks" that formed 
this halo. The simulated Mv-[Fe/H] relation of MPCs can be 
compared with the observed "blue tilt". The least-square fit to 
the five simulation data shows Z oc L -2 01 . 



elliptical galaxies (e.g., Djorgovski & Davis 1987; Capelato 
et al. 1995). R c and a of GCSs are significantly different 
from the half-mass radii and central velocity dispersions of 
dark matter halos that host the GCSs. In the present sim- 
ulations, MPCs within galaxy-scale halos at z — trace 
the particles initially within the central region of subhalos 
virialized before z — z t r Un . On the other hand, dark matter 
halo particles at z — trace all the particle from subhalos 
virialized at every z. Therefore, R e and a of GCSs are quite 
different from the half-mass radii and central dispersions of 
dark matter halos that follow the virial relation. 

Fig. 10 shows the best fit to the three properties of R e , 
h, and a for GCSs (i.e., R c = C*pp/ C " ' 46 cr ' 60 ). By assum- 
ing that R c oc 7 C QFP x ct' 3fp , the best values of a^p and /3fp 
are chosen such that the dispersion along the assumed line 
(or plane) with app and (3fp become the smallest. The de- 
rived relation of R c = CppI c ~ 0A6 a 0S0 is quite different from 
the virial relation of R c = C*vp/c~V 2 , which suggests struc- 
tural and kinematical non-homology in the simulated GCSs. 
The more significant deviation of /3fp from the virial relation 
suggests that the estimated line-of-sight velocity dispersions 
of GCSs can be quite different from (or significant smaller 
than) the central velocity dispersions of dark matter halos 
required for virial equilibrium. Since these three properties 
of GCSs are feasible to derive observationally, it is an in- 
teresting future study to compare these simulated scaling 
relations with the corresponding observational ones. 
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Figure 12. The same as Fig. 11 but for / gn = 0.0001. The 
least-square fit to the five simulation data shows Z oc L~ 2 ' 69 . 
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Figure 13. The number of galactic nuclei vs halo mass for 
halos with Ztrun = 10. 



3.3 Blue Tilts 

Next we describe the addition of accreted stellar galactic 
nuclei (GN) to our simulated GCSs. Fig. 11 shows that the 
GCS in a massive halo of M h = 3.0 x 10 13 M Q for the model 
in which the fraction of the stellar galaxy mass in nuclei is 
/ gn = 0.001. The figure shows a clear [Fe/H]-M v relation 
in the more luminous MPCs (Mv < —9 mag). The essen- 
tial reason for the simulated blue tilt is that this halo was 
formed from a large number of subhalos (N h — 62) that 
were virialized before z t r Un and thus had GN that followed 
a mass-metallicity relation. Whether or not the blue tilt in 
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a halo at z = can be clearly seen is determined by how 
many subhalos with GN particles merged to form the halo 
in the present simulations. A least square fits to the simu- 
lation data show that M v = -10.25 - 1.24 x [Fe/H] (i.e., 
Z oc L 2 ' 01 ), which means that the simulated blue tilt is sig- 
nificantly steeper than the observed one (e.g., Z oc L ' 55 ; 
Harris et al. 2006). 

Fig. 12 shows that if / gn = 0.0001, which is an order 
of magnitude smaller than the model in Fig. 11, the blue 
tilt is much less clearly seen, because low-luminosity GN 
particles can not be distinguished from GC particles in the 
GCS. The results in Figs. 11 and 12 imply that for a blue 
tilt to be clearly seen at z = 0, GN at z = 10 should be 
as massive as ~ 0.1 % of their stellar components in their 
low-mass hosts. We confirm that if / gn = 0.005, a blue tilt 
can be clearly seen but there are too many very luminous 
(My < — 12 mag) MPCs in galaxy-scale halos. The least 
square fits to the simulation data show that Mv = —8.22 — 

0. 93 x [Fe/H] (i.e., Z oc L 2 ' 69 ), which is steeper than the blue 
tilt in the model with / n = 0.001 (and thus less consistent 
with observations). These results suggest that the presence 
of blue tilts gives some constraints on the destruction of 
low-mass nucleated galaxies around galaxies between z trU n 
and z = 0. The dependences of the slopes of the simulated 
blue tilts on the initial GC luminosity function are given in 
Appendix E. 

Fig. 13 shows that more massive halos will have a larger 
number of GN particles within their GCSs and thus massive 
galaxies at z = are more likely to have a contribution 
from disrupted nuclei to their observed blue tilts. The sig- 
nificantly smaller dispersion in JV gn for more massive halos 
with M h > 10 13 M Q in Fig. 13 suggests that this effect is 
most pronounced for GCSs in the most luminous galaxies, 

1. e. those at the centers of groups and clusters. 

The overall distribution of MPCs in the M v -[Fe/H] 
plane does not reproduce well the observations for GCSs 
with blue tilts (Harris et al. 2006; Strader et al. 2006; Spitler 
ct al. 2006). The stripped nuclei, in this prescription, arc 
therefore unlikely to be the sole origin of the observed blue 
tilts (although they may still make a small contribution). 



4 DISCUSSION 

4.1 The origin of the L — Z sc relation 

The present numerical simulations are the first to demon- 
strate that an initial steep L — Z gc relation at high z can 
become significantly flatter due to hierarchical merging of 
galaxies with GCSs. The slope of the final L — Z gc rela- 
tion GCSs at z = is similar to the observed one (Strader 
et al. 2004; P06) in models with 2trun = 10 and constant 
M/L at Ztrun- Cote et al. (1998) investigated correlations be- 
tween metal-poor GCSs and their host galaxy luminosities 
(L) in their models of GCS formation via merging/accretion 
of MPCs from dwarfs with a power-law galaxy luminos- 
ity function (LF) of slope a. They showed that there is a 
weaker correlation between the MPCs and L in the mod- 
els with a steep LF (a ~ —1.8 in their Figure 7). Recent 
high-resolution cosmological simulations on the mass func- 
tion of low-mass halos have demonstrated that the slope of 
the mass function is significantly steeper at higher redshifts 



(e.g., Yahagi et al. 2004). We suggest that the steeper mass 
function of halos at high redshifts can also be an important 
factor for the origin of the observed slope of the L — Z gc 
relation. 

The present study assumed that the initial L — Z gc rela- 
tion at z = ztrun is similar to the L — Zs relation observed in 
Local Group dSphs today (e.g., Dekel & Silk 1984; Lotz et al. 
2004). Such a steep L — Z gc relation at z = .ztrun is important 
for the present model, because the initial L — Z gc relation 
becomes significantly flatter between z = Zt run and z — 0. 
It is, however, observationally and theoretically unclear how 
low-mass galaxies at high z achieve such a steep relation 
and what physical mechanisms are responsible. Kravtsov & 
Gnedin (2005) found that galaxies in their cosmological sim- 
ulations exhibit a strong correlation between the stellar mass 
(Ms) and the average metallicity (Zs) of stars (described as 
Zs oc Ms ' 5 ) similar to the observed one in dwarf galaxies 
(Dekel & Woo 2003). Their results imply that GCSs in low- 
mass galaxies at high z could also have a steep Z gc — Ms 
relation similar to that of the field stars. 

It should be stressed here that the models with an ini- 
tially steep L — Z gc relation at z = Ztrun better reproduce 
observations only for constant M/L at z = z tIun : models 
with an initially flatter L — Z gc relation can also better ex- 
plain observations for variable M/L at z = z tIun , though the 
required L — Z gc relation was not predicted by previous the- 
oretical studies. Owing to lack of observational data sets for 
the L — Z gc relation at z = 2trun, it is currently difficult to 
make a robust conclusion as to whether the constant M/L 
models with a steep L — Z gc relation are better. 

4.2 The origin of the blue tilt 

Using the ACS on-board HST, Strader et al. (2006) and 
Harris et al. (2006) found that luminous blue GCs (i.e. 
MPCs) reveal a trend of having redder colors in giant el- 
lipticals. This "blue tilt" was interpreted as MPCs having 
a metallicity-luminosity relation of Z oc L ' 55 . A possible 
exception was NGC 4472. Spitler et al. (2006) subsequently 
showed that this trend is also true in the Sombrero spi- 
ral galaxy and may extend to lower GC luminosities with 
a somewhat shallower slope than derived by Harris et al. 
(2006) and Strader et al. (2006). In each of these ACS stud- 
ies the MRCs showed no corresponding trend. These obser- 
vations require that any theoretical models of GC formation 
should explain the origin of the metallicity-luminosity rela- 
tion for individual GCs, and also the apparent absence of a 
similar relation for MRCs (e.g., Mieske et al. 2006). 

Here we have demonstrated that luminous MPCs in 
more massive galaxies show metallicity-luminosity relations 
in that brighter GCs have redder colors — a trend qualita- 
tively similar to the observed blue tilts. In our simulations 
the reason for the "simulated blue tilts" is that luminous 
MPCs originate from stellar galactic nuclei (GN) of the more 
massive nucleated galaxies with a luminosity-metallicity re- 
lation. The present study thus suggests that galaxies which 
experienced a larger number of accretion/merging events of 
nucleated low-mass galaxies are more likely to show a blue 
tilt. Since the frequency of accretion/merging events is gen- 
erally higher in more massive galaxies, the blue tilts are 
more likely to be present in high luminosity than low lumi- 
nosity galaxies. Thus not all galaxies are expected to have 



Formation of GCSs 11 



blue tilts. Since these simulated GCSs are composed both of 
stripped normal GCs (with no initial mass-metallicity rela- 
tion) and of stripped GN with a mass-metallicity relation, 
their detection in observational data may depend on their 
relative numbers. 

Although we find clear blue tilts in our simulations, it 
is not clear what role stripped galactic nuclei play in the ob- 
served blue tilts to date. For example, Spitler et al. (2006) 
find the blue tilt is not simply a feature of luminous MPCs 
but seems to extend down to luminosity function as far as 
is measurable. Furthermore the MPC luminosity function is 
well fit by a standard Gaussian, or t 5 function, even at the 
luminous end. Thus the tilt does not appear to be due to an 
additional population but rather a 'tilting' of the existing 
MPC system. More deep data on a variety of systems are 
required to confirm these trends. A direct comparison be- 
tween simulations and observational data needs to account 
for many factors, such as the number of GN vs bona fide 
GCs, contamination rates, the presence of dust and measure- 
ment errors (which tend to cause spreads in both magnitude 
and colour). 

The observed colors and metallicities of GN in dwarfs 
(Lotz et al. 2004; Mieske et al 2006) appear to suggest most 
GN can not be metal-rich clusters as they have lower metal- 
licities. The apparent absence of a MRC color-magnitude 
relation may therefore reflect the fact that most MRCs do 
not originate from GN. Our simulated blue tilt, due to the 
stripped nuclei of (defunct) low-mass galaxies, is one of 
many possible scenarios for the origin of the blue tilt. We 
plan to discuss the origin of the blue tilt in a wider context 
of galaxy and GC formation in a future paper. 

4.3 The Fundamental Plane of metal-poor GCSs 

The scaling relations between properties of elliptical galax- 
ies, such as "the Fundamental Plane" (FP), have long been 
suggested provide valuable information on formation and 
evolution histories of elliptical galaxies (e.g., Djorgovski & 
Davis 1987). Theoretical studies based mostly on numerical 
simulations suggested that structural and kinematic non- 
homology inferred from the FP are closely associated with 
dynamics of galaxy merging, star formation histories depen- 
dent on galaxy masses, and collapse dynamics in elliptical 
galaxy formation (Capelato et al. 1995; Bekki 1998; Dantas 
et al. 2003). The present study has shown that the scal- 
ing relation (R e = CFp/ e _0 ' 46 o" ' 60 ) between R c (projected 
half- number radii), I c (surface number density at R e ), and 
a (velocity dispersion) of MPCs in our models is signifi- 
cantly different from that expected from the virial theorem 
(R c — Cvp/c^cr 2 ). Although this is due partly to the way 
we estimate a for a GCS in our models, the significant dif- 
ference between the exponents of the I e terms between the 
virial and the simulated relations strongly suggests that the 
deviation has a physical meaning. 

MPCs are assumed to form in the central regions of 
low-mass halos that are virialized before z = ztrun- The for- 
mation of MPCs is biased toward high density peaks of the 
primordial matter distribution in the adopted CDM model. 
As a natural result of this, MPCs in galaxies at z = show 
structural and kinematical properties different from those 
of the background dark matter halos, the scaling relation 
of which should be similar to the virial relation. Recent nu- 



merical simulations with truncation of GC formation have 
also demonstrated that the dynamical properties of MPCs 
can be significantly different from those of dark matter ha- 
los (Santos 2003; Bekki 2005; Moore et al. 2006). Although 
physical explanations for the origin of the slope of the FP 
in MPCs are yet to be provided, we suggest that the pres- 
ence of the FP for MPCs can be one characteristic of biased 
formation of MPCs at high z. 

A growing number of spectroscopic observations have 
revealed kinematic properties of GCSs, such as radial profiles 
of velocity dispersions and rotational velocities, maximum 
V/a, and kinematical differences between MPCs and MRCs, 
in early- type galaxies (e.g., Kissler-Patig & Gebhardt 1998; 
Cohen 2000; Zepf et al. 2001; Cote et al. 2001, 2003; Beasley 
et al. 2004; Peng et al. 2004; Richtler et al. 2004; Pierce et 
al. 2006; Romanowsky 2006). Recent photometric studies of 
GCSs by the HST and ground-based telescopes have also 
revealed structural properties of GCSs in galaxies with dif- 
ferent Hubble types (e.g., Rhode & Zepf 2004; P06; Spitler 
et al. 2006). Accordingly it is worthwhile for future obser- 
vational studies to investigate the FP of MPC systems. It is 
also an interesting observational study to investigate corre- 
lations between R c (or 7 C ) of MPCs and the total luminosity 
of their host galaxy, because such correlations are more fea- 
sible to derive observationally and thus can be compared 
more readily with the predicted ones. 



5 CONCLUSIONS 

We have investigated the physical properties of metal-poor 
GCSs in galaxies using high-resolution cosmological simula- 
tions based on a ACDM model. We particularly investigated 
the L - Z gc relation, the My - [Fe/H] relation, and MPC 
scaling relations in galaxies with different masses. We as- 
sumed that (i) MPC formation is truncated at z = z tIun , 

(ii) galaxies at z = z tIun contain MPCs and galactic nu- 
clei (GN) with masses proportional to their host galaxies, 

(iii) galaxies initially have a steep relation at z — 2 t run, (iv) 
galaxies have a certain M/L at z = z trU n, and (v) MPCs 
have a luminosity function the same as that observed today. 
We also assumed that both MPCs and GNs can be identi- 
fied as MPCs in galaxies at z = 0, if they are within virial 
radii of the dark matter halos. We summarize our principle 
results of the models as follows. 

(1) The original L — Z gc at z trU n becomes significantly 
flatter by z — due to hierarchical merging of lower-mass 
galaxies. The original Z gc oc L 5 at z = 10 in the model 
with a constant M/L at Ztrun becomes Z gc oc L 015 at z = 0, 
which is consistent with the latest observational results. The 
origin of the flattening of the L — Z gc relation reflects the fact 
that Z gc of GCSs in more massive galaxies at z = is lower 
than the progenitor halos at z trU n due to a large number of 
mergers and accretion of low-mass halos of GCSs with lower 
Z gc . A flatter L — Z gc relation at z — 10 {Z gc oc L ' 25 ) is 
required to explain the observed relation at z — in the 
models with a variable M/L. 

(2) The final L — Z gc relation at z = depends on z tIU n 
such that it is steeper for lower z tm n. Models with constant 
M/L ratios (i.e., M h /L = 10) at z — -^trun can better explain 
the observed L — Z gc relation at z — for a given z = z trun 
and the adopted initial L — Z gc relation of Z gc oc L ' 5 at 



12 K. Bekki, H. Yahagi, and D. A. Forbes 



z = Ztrun , though most models show a flattening of the orig- 
inal L — Zg C relation. The essential reason for the positive 
correlation between L and Z gc at z = is that the more mas- 
sive galaxies at z — are formed from hierarchical merging 
of a larger number of more massive building blocks contain- 
ing GCs with higher metallicities. Thus merging does not 
completely erase the original L — Zg C relation. 

(3) MPC systems in these models clearly show scaling 
relations between galaxy mass M s (or halo mass Mh), -Re, 
J e , and g. We find R e oc M s °- 85 , 7 e ,ge oc Ms" - 97 , a oc M s °' 48 
in the models with z tIU n = 10. Thus we find a "Fundamental 
Plane" for MPCs, which can be described as Re oc 7 e ~°' 5 f °' 6 
in the models with zt lun = 10. We note that such scalings 
are not consistent with the virial relation of R e oc I c _1 <t 2 . 
The R c — M s relation depends on Zt run , which suggests that 
the scaling relation can give some constraints on Ztrun- We 
also find that Re oc M h ' 57 and a oc M h ' 32 . 

(4) Luminous MPCs show a correlation between My 
and [Fe/H] if these MPCs originate from nuclei of low-mass 
galaxies at high z. The correlation, which is referred to as a 
"blue tilt" , can be more clearly seen in more massive galax- 
ies with Mh ~ 1O 13 M0. This is mainly because more mas- 
sive galaxies are formed from a larger number of nucleated 
galaxies virialized by z = z tr un- Observations of blue tilts 
only loosely resemble our simulated ones suggesting stripped 
nuclei make a small contribution at most. 
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APPENDIX A: DEPENDENCES OF M H Zgc 
RELATIONS ON Ztrun. 

Fig. Al shows that the initially steep Mh — Z gc relation at 
z — 6 becomes flatter in the model with z trull = 6, though 
the degree of the flattening is less significant in this model 
than in the model with ztrun = 10. This reflects the fact that 
galaxies experienced a smaller number of galaxy merging 
(which can flatten the Mh — Z gc relation) between z — z trun 
and z — in this model. Fig. A2 shows that the Mh — 
Z gc relation is almost flat in the model with ztrun = 15, 
which means that hierarchical merging of low-mass galaxies 
between ztrun and z — almost completely erases out the 
original steep Mh — Z gc relation: this simulated Mh — Z gc at 
z = is inconsistent with the observed one. These results 
shown in Figs. Al and A2 thus suggest that the degree of the 
flattening in the Mh — Z gc relation depends on Ztrun in the 
sense that models with earlier truncation of GC formation 
show the larger degree of the flattening. The absence of low- 
mass galaxies with GCs for log 10 Mh > 11.4 (M@) in the 
model with Zt rU n = 15 is inconsistent with observations. 



APPENDIX B: DEPENDENCES OF M s - Z GC 
RELATIONS ON M 3 /L. 

Fig. Bl shows that the original L — Z gc relation becomes sig- 
nificantly flatter between ztrun and z = in the model with 
M s /L = 5. The locations of galaxies on the M s -Z gc plane 
are shifted to the righter directions owing to the larger stel- 
lar masses of galaxies in this model with higher M s /L. As 
a result of this, the simulated Z gc for a given M s is slightly 
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Figure Al. The same as Figure 2 but for the model with zt run = 
6. 
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Figure A2. The same as Figure 2 but for the model with zt run = 
15. 

smaller than the observed one. The change in M&/L how- 
ever does not change the final slope in the L — Z sc relation in 
comparison with the model with M s /L = 1. The flattening 
of the L — Z gc relation between 2 tnm and z = does not 
depend on M s /L, as long as we adopt assumptions of (i) a 
constant M s /L at z — z trun and (ii) an initial L — Z gc rela- 
tion of Z gc oc L 0,5 . This means that if we adopt a slightly 
higher value of /3 gc for the higher M s /L model, the observed 
relation can be better reproduced under the above assump- 
tions. Thus high values of /3 gc are required in the models 
with higher M s /L and with the above assumptions for the 
observation to be well reproduced. 
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Figure Bl. The same as Figure 3 but for the model with Figure CI. The same as Figure 3 but for the model with ztrun = 

Ms/L = 5. 6. 



APPENDIX C: DEPENDENCES OF Ms - Z GC 
RELATIONS ON Ztrun. 

Figs. CI and C2 show the simulated M s — Z gc relations (thus 
L — Z gc relations for the fixed M s /L ratios) for the models 
with 6 and 15, respectively. It is clear from these 

figures and Fig. 5 that the model with z trun = 10 can better 
reproduce the observed flat L — Z gc relation (or M a — Z gc 
one). The model with ztrun = 15 shows the L — Z sc rela- 
tion flatter than the observed one whereas the model with 
ztrun = 6 shows the steeper one. In the model with ztrun = 15 
(ztrun = 6), a larger (smaller) number of merger events for 
a longer (shorter) time scale between z — Ztrun and z = 
can flatten the original L — Z gc relation at z — Ztrun to a 
larger (smaller) extent. This is the essential reason why the 
simulated L — Z gc relation depends on ztrun. 

Thus the growth of galaxies via hierarchical merg- 
ing/accretion of low-mass halos with and without MPCs be- 
tween z — ztrun and z — can flatten the original L — Z gc at 
2 = £trun owing to merging/accretion of MPCs with different 
metallicities. If Zt tun corresponds to the epoch of reioniza- 
tion (zreion) and it is determined from future observations, 
then the results of numerical simulations shown in Figs. 6, 
CI, and C2 imply that we can infer the original L — Z sc at 
z — ^rcion by comparing the simulations with the observed 
L — Zg C at z — 0. 



APPENDIX D: DEPENDENCES OF SCALING 
RELATIONS ON Ztrun- 

Fig. Dl shows that the dynamical correlations are different 
between different Ztrun, in particular, for M s -dependences 
of I c . the M s -dependence of R e is steeper for higher ztrun 
whereas the M s -dependence of a is not so different between 
different ztrun - The slope in the M s -dependence of I c is nega- 
tive for ztrun = 10 and 15 and positive for ztrun = 6. The de- 
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Figure C2. The same as Figure 3 but for the model with ztrun = 
15. 



rived differences in the M s -dependences of dynamical prop- 
erties of GCSs between different z trun suggest that future 
observational studies on the Ms-dependences can provide 
some constraints on ztrun- 

Structural properties of GCSs have been suggested to 
provide constraints on ztrun (Santos 2003; Bekki 2005; Moore 
et al. 2006). The present study suggested that if (i) the 
adopted L — Z gc relation at z = Ztrun {Z gc oc L ' 5 ) is reason- 
able and realistic and (ii) M/L does not depend on galaxy 
masses at z = ztrun, the observed L — Z gc relation at z — 
could also provide constraints on z tmn . It is therefore an 
important observational test whether Ztrun derived from the 
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Figure Dl. Dependences of linc-of-sight velocity dispersions 
(er, top), effective surface number densities (J e , middle), and half- 
number radii (R c , bottom) of MPCs on M B at z = for the 
model with ztrun = 6 (solid), ztrun = 10 (dotted), and ztrun = 15 
(dashed). Each of these lines represent the least-square fit to the 
simulation data for each of the three Ms-dependences in each 
model. 



observed M s -dependences is consistent with z trU n derived 
from the structural and chemical properties of GCSs. 
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Figure El. The same as Fig. 11 but for the model in which 
maximum GC masses (m max ) in halos are assumed to correlate 
with their host halo masses. 



overall distribution of MPCs on the Mv-[Fe/H] in this model 
appears to be only slightly more similar to the observed blue 
tilt than the model shown in Fig. 11. The least square fits to 
the simulation data show that M v = -10.74-1.62 x [Fe/H] 
(i.e., Z oc L 1 ' 55 ). This Z oc L 1 ' 55 relation is however still sig- 
nificantly steeper than the observed one of Z oc L ' 55 (Harris 
ct al. 2006). Although the models with steeper m max — Mh 
relations can show flatter Z — L relations, the simulated 
dispersions in My for given metallicity bins (in particular, 
for [Fe/H] > —1.5) are too large to be consistent with ob- 
servations. These suggest that the models with m max — Mh 
relations still miss some important ingredients of GC forma- 
tion. We discuss these problems in our forthcoming papers 
(Bckki ct al. 2007, in preparation). 



APPENDIX E: THE INFLUENCES OF 
MAXIMUM GC MASSES ON THE BLUE TILT 

We adopted an assumption of a universal GC mass (lumi- 
nosity) function with the same lower (mi ow ) and upper (or 
luminosity) mass cut off (m upp ) for all halos in the present 
study. Since recent observations have suggested a possible 
maximum mass m max of star clusters in galaxies (e.g., Gieles 
ct al. 2006) and correlations between m max and physical 
properties of their host galaxies, such as star formation rates 
(e.g., Larsen & Richtler 2000), the above assumption on the 
fixed m up p could be less realistic. We thus investigated the 
models with m max — Mh relations (i.e., m upp — Mh) suggested 
by numerical simulations of Kravtsov & Gnedin 2005. 

Fig. El shows the result of the model with M h = 3.0 x 
10 13 M Q and N gn = 62 for m max = 2.9x 1O 6 M ( T5 £^) 1 ' 29 . 
Owing to the introduction of the m max — Mh relation, the 



